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Abstract 


In wind-diesel stand-alone power systems, the disturbances like random nature of wind power, turbulent wind, sudden changes in load demand 
and the wind park disconnection effect continuously the system voltage and frequency. The satisfactory operation of such a system is not an easy 
task and the control design has to take in to account all these subtleties. For maintaining the power quality, generally, a short-term energy storage 
device is used. In this paper, the performance of a wind-diesel system associated with a superconducting magnetic energy storage (SMES) system 
is studied. The effect of installing SMES at wind park bus/load bus, on the system performance is investigated. To control the exchange of real and 
reactive powers between the SMES unit and the wind-diesel system, a control strategy based on fuzzy logic is proposed. The dynamic models of 
the hybrid power system for most common scenarios are developed and the results presented. 


€ 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


In most of the countries worldwide, there are remote com- 
munities, which are not served by the main electrical grids and 
may never be connected to them because of economical reasons. 
Diesel generators are mainly used to serve them. But numerous 
remote areas have significant wind energy potential. Wind-diesel 
hybrid power systems are increasingly being seen as a viable 
and economical alternative to grid supply for such locations. 
One major advantage of a wind energy based stand-alone power 
system is that it reduces diesel fuel consumption. However, such 
stand-alone power systems are characterized by low inertia and 
poor reactive power support and therefore face wide voltage and 
frequency variations due to changes in load demand. The inte- 
gration of wind turbines in such systems provokes additional 
power quality problems in the form of frequency and voltage 
fluctuations, due to rapidly changing nature of the wind [1-3]. 

The existing method used to solve these problems involves 
installation of batteries. The batteries are however constrained 
by their load cycles. There are some environmental concerns 
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related to battery storage due to toxic gas generation during 
battery charge/discharge operation. The disposal of hazardous 
materials presents some disposal problems. On the other hand, 
the technical improvements in SMES and a better knowledge of 
dealing with cryogenic systems have made the SMES units to 
penetrate the market very fast [4-7]. The dynamic performance 
of an SMES system is far superior to other storage technolo- 
gies. The shorter response time is the leading advantage and 
at the same time, system operation and life time are not influ- 
enced by the number of duty cycles or the depth of discharge as 
is the case for classical batteries. The estimated life of a typical 
SMES system is about 30 years. With the on going fast evolution 
in power electronics and SMES technology the estimated cost 
within 10 years is expected to be almost 25% of present cost. 
SMES systems are compact, self-contained and highly mobile. 
They contain no hazardous chemicals, produce no flammable 
gases and can be kept at remote locations. The SMES unit can 
be controlled to produce both active and reactive power mod- 
ulations in the four quadrants of the P-Q plane. It can thus 
be concluded that the SMES system therefore perfectly fits the 
wind-diesel system for power quality improvement. 

To exchange the real and reactive powers between the SMES 
unit and an ac bus of the wind-diesel power system suitable 
control strategies are required. The parameters of a wind-diesel 
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Nomenclature 


Eq, Eq’ synchronous generator voltage behind subtran- 


Eb, E, 


Eta 


Q 


sient reactances 

real and imaginary components of induction 
machine voltage behind transient reactance(p.u.) 
synchronous generator field voltage (p.u.) 
synchronous generator voltage behind transient 
reactance (p.u.) 

frequency(Hz) at load bus/wind park bus 

total wind turbine and wind generator inertia 
constant (s) 

direct and quadrature axis components of 
synchronous machine injected current 

direct and quadrature axis components of 
induction machine injected current 


IDsm/Qsm D-Q components of SMES current 


Ism 


current supplied by SMES 

nominal current of SMES coil 

SMES coil inductance 

real power demanded by SMES unit 

real power supplied by SMES unit 

real power inputted to SMES unit 

reactive power demanded by SMES unit 
reactive power supplied by SMES unit 
reactive power inputted to SMES unit 

stator resistance of synchronous machine (p.u.) 
stator resistance of induction machine (p.u.) 
induction machine slip 

electromagnetic torque developed by induction 
machine (p.u.) 

mechanical torque developed by wind turbine 
(p.u.) 

converter time delay 

Electromagnetic torque developed by 
synchronous machine (p.u.) 

mechanical torque developed by diesel engine 
(p.u.) 

Induction generator rotor open circuit time 
constant (s) 

direct and quadrature axis components of 
synchronous machine terminal voltage (p.u.) 
real and imaginary components of induction 
machine terminal voltage (p.u.) 

SMES coil voltage 

dc voltage of each converter with firing anglr 
equal to zero 

magnitude of voltage of induction generator 
bus/load bus (p.u.) 

induction generator open circuit reactance (p.u.) 
induction generator open circuit transient 
reactance (p.u.) 

synchronous generator reactance of direct and 
quadrature axis(p.u.) 


X3, X ü synchronous generator transient reactance of 
direct and quadrature axis (p.u.) 

X'i, XG synchronous generator subtransient reactance of 
direct and quadrature axis (p.u.) 


Greek letters 
OA wind generator angular velocity (p.u.) 
wo synchronous angular velocity 


power system change with operating conditions. Thus, the best 
performance cannot be achieved by fixed gain controllers. To 
improve the performance of a wind-diesel power system for 
different operating conditions, self-tuning controllers were pro- 
posed [8]. A major disadvantage of the self-tuning control 
scheme is the need to identify system model parameters in real 
time, which is very time consuming. In order to get rid of the 
computational burden associated with self-tuning control, a dif- 
ferent control strategy based on fuzzy logic is proposed, in this 
paper. The advantage with the fuzzy logic controllers is their 
ability to effect non-linear control actions, which are desirable 
for better system dynamic performance. There is no need for 
parameter identification and these controllers are less depen- 
denton rigorous mathematical models. The proposed fuzzy logic 
control scheme is very simple and has a small number of rules. 
The input signals for this scheme are SMES coil current, fre- 
quency and voltage deviations at wind park bus/load bus. Various 
models are developed for various perturbations and system con- 
figurations. The effect of associating the SMES unit with the 
wind-diesel system is then investigated using simulation stud- 
ies. In one of the configurations of the hybrid system, the SMES 
unit is connected to load bus and the performance of this con- 
figuration is compared with that of another case in which the 
SMES unit is connected to wind park bus. 


2. System under study 


Fig. 1 shows the system under study [9], which consists of two 
synchronous machines driven by diesel engines, ten induction 
machines driven by stall regulated wind turbines, transmission 
network, capacitor banks and a load. In one of the system config- 
urations, the proposed SMES is connected to the wind park bus 
while in another configuration it interacts with the wind-diesel 
system through the load bus. 


3. Description of system components 
3.1. Wind turbine 


Each stall regulated wind turbine and its coupled induction 
generator is represented by a lumped single mass, so that equa- 
tion of motion can be written as: 


d = : T, T, 1 
a On) = Ha Am — TAe) (1) 
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Fig. 1. Wind-diesel-SMES System. 


3.2. Induction generators 


Each induction machine coupled to a wind turbine is an asyn- 
chronous machine operating with a negative slip and electro- 
magnetic torque. The model of each machine in synchronously 
rotating (D-Q) reference frame is expressed as: 

Stator equations (in p.u.): 


Vp = Ep +r3lp — X'IQ (2) 

Vo = Eg  rsslo + X Ip (3) 
Rotor equations (in p.u.): 

d (gi j= |e — (X — X’ lo] + Sao E; (4) 

ap ee oe 

d / 1 u / 

gEV = 7! — (X — X Ip] + Sao Ep (5) 


Electromagnetic torque equation: 


TAe — EpIp + Eolo (6) 
3.3. Capacitor compensation 


For the induction generator, driven by a wind turbine, it is 
a well-known fact that it can deliver only active power and 
requires reactive power from the grid. The requirement can be 
quite significant even at no-load. The induction generator has 
a reactive power requirement, which varies with active power 
generated by it. In this paper, the 7596 of the reactive power 
requirement is compensated with the shunt capacitors and the 
remaining portion of the reactive power requirement is drawn 
from the network. The capacitor banks installed for compen- 
sation of reactive power demanded by induction machines are 


modeled by incorporating their equivalent reactance value in the 
network admittance matrix. 


3.4. Synchronous machines and their associated devices 


Each synchronous machine driven by a diesel engine unit is 
represented by voltage and torque equations in machine (d-q) 
reference frame as: 

Stator equations (in p.u.): 


Va = Ej = si là = X, Iq (7) 


Va = Ej — rsilq — Xola (8) 


(i—1,2) 
Rotor equations (in p.u.): 


MD D= E [Eq(Xq — Xp] (9) 
d 1 Xa- X1 Xa — X, 

T(E) ee |En- —— Eq ug (10) 
dt - Tio X-X 9 X4—-Xj 4 

d 

gE = = A - E+ (Xa - XH] (11) 


Electromagnetic torque equation is: 


Tp. = Ejla + Ela -(X4— Xalq (12) 

The standard models are used for the diesel engines and their 
associated speed governors [8,9]. The AVRs used are of IEEE 
typel [8,9]. 
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Fig. 2. SMES unit dynamics and control. 


3.5. Transmission network 


With the transmission line characteristics (Rr, XT) per unit 
length known, they could be multiplied by the line length [8,9] 
and the resulting values considered concentrated over transmis- 
sion line. 


3.6. Load 


In a particular configuration, the connected load is repre- 
sented by a constant admittance and its admittance value is 
incorporated in the network admittance matrix. 


3.7. SMES system 


The various forcing functions (perturbations) in a wind-diesel 
system upset the real and reactive power balance and thus affect 
the system frequency and voltage. Before the diesel engine gov- 
ernors and AVRs will act, a control device like SMES which 
is capable of modulating fast, the active and reactive powers is 
an appropriate candidate for power quality improvement of the 
system. When a superconducting coil is connected through a 


GTO/IGBT based, 12-pulse converter to a three-phase ac bus, 
it can modulate real and reactive powers. The control of the fir- 
ing angles a; and oo of the cascaded six-pulse converters gives 
the SMES unit capability to exchange (inject/absorb) active and 
reactive powers, independently with the system [7]. For initial 
charging of the SMES coil the converter voltage is held pos- 
itive at a suitable constant value. When the inductor current 
reaches its nominal value /9,. the SMES unit is connected to 
wind park/load bus and left floating by adjusting o to 90°, a2 
to 270°. Under these conditions, there is no exchange of real and 
reactive powers between the SMES unit and the wind-diesel sys- 
tem. As soon as the controllers associated with the SMES unit 
experience a change in voltage and frequency at the bus SMES 
unit is connected to the wind-diesel system, the firing angles 
of the converters are controlled to modulate active and reactive 
powers. The block diagram of Fig. 2 depicts the modeling and 
the control of SMES unit. Details about this block diagram are 
given in the following sub sections. 


3.7.1. Fuzzy logic controllers for SMES 
Two fuzzy logic controllers, P- and Q-controller generate, 
Pgem and Qdem, which represent, respectively, the active and 
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Fig. 3. Membership functions for P and Q controllers. 


reactive powers required to be supplied from the SMES system. 
The P-controller uses a mixture of ac bus frequency deviation 
and SMES coil current deviation signals. For continuous control, 
it is necessary that SMES coil current has a tendency to attain 
its nominal value, that is why an appropriate portion of this sig- 
nal is added to frequency deviation signal. The Q-controller on 
the other hand requires the ac bus voltage deviation for gener- 
ating Qaem. Triangular type of membership functions are used 
for input as well as output fuzzy sets as shown in Fig. 3. For 
the P-controller, the universe for discourse for the input (X) 
units and output Pgem are chosen as (—0.3 units, 40.3 units) and 
(—30 kW, +30 kW), respectively. For Q-controller, the universe 
of discourse for input (Y) units and the output Qdem are cho- 
sen as (—0.02 p.u. and +0.02 p.u) and (30 KVAR, +30 KVAR), 
respectively. The fuzzy rules are determined from Table 1. Seven 
fuzzy sets are considered for input variables and five sets for out- 
put variables and the height defuzzification method is used for 
computing control signals. 


3.7.2. Calculation of firing angles 
The firing angles of two cascaded six-pulse converters under 


four-quadrant operation are calculated from Toon (2—Pgem) and 


Jer (=—Qdem) as per the following equations [7]: 
P! 
aj = cos™! ( — m. 
dem dem 
/2 D 
=j P dem + Q dem 
+ cos — aro (13) 
2V2 Tem 
P! 
az = cos™! ( — m 
dem dem 
/2 12 
zj P dem + Q dem 
— cos"! (ALLTEL <m (14) 
ay Ism 
Table 1 
Rule base for active/reactive power generators 
XorY LN MN SN Z SP MP LP 
Paem OF Qdem LP LP SP Z SN LN LN 


3.7.3. SMES unit dynamics 

The actual values of active and reactive powers, Psm and Qsm, 
injected by SMES at the ac bus are calculated after taking into 
account the converter delay and the SMES dynamics presented 
in Fig. 2. In this figure, the outputs of the blocks representing the 
converter delay are identified as @ and à?. The various equations 
governing the SMES unit are then. 

Voltage across the SMES coil 


Vim = vo (cos à, + cos à2) (15) 
Current through the SMES coil 


Ism = 


t 
I Vem dt + 7, (16) 
ti 


sm 0 


Active power injected by SMES unit at ac bus 

Pm = — VÀ, Is (cos & + cos 2) (17) 
Reactive power injected by SMES at ac bus 

Osm = — VO, Ia (siná; + sin &2) (18) 


Once these equations are cast in the block diagram form, the 
middle portion of Fig. 2 (labeled as SMES unit dynamics) is 
obtained. 


3.7.4. SMES unit model in D-Q rotating frame 

The D and Q components of the current injected by SMES 
unit at the ac bus, in response to calculated firing angles a; and 
a2, are calculated as per the following equations: 


Vp P, Vo Qs 

IDsm = on +4 Ofen (19) 
V5 + Vo V5 + Vo 
Vo Psm Vp Osm 

1Qsm = 2 E 2 2 : 2 Q0) 
V5 + Vo V5 + Vo 


In Eqs. (19) and (20), Vp and Vg are D-Q components of ac bus 
voltage, where SMES is connected. The complete SMES unit 
dynamics and control scheme is given in Fig. 2. In this figure, Ke 
represents the SMES coil current feedback gain the importance 
of which is already given in Section 3.7.1. 


4. Compact network-machine model 

In general, the network equations for the hybrid wind-diesel 
SMES system can be written as: 
I = [YBuslV Q1) 


where [Ypus] is the admittance matrix of the system, formulated 
with equivalent admittances of the capacitor banks, transform- 
ers, and loads incorporated. In the partitioned form, the Eq. (21) 
can be written as: 


In) _ | na] Ward | | Va (22) 
0 Yml p he 
where subscript n denotes the total number of network nodes at 
which (synchronous and induction) generators and SMES unit 
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are connected. The subscript r on the other hand denotes the 
remaining nodes. The above equation can be simplified to take 
the form of (23) 


Vn = [ZTT] In (23) 
where 
[Yrr] = ([Ya4] — [Ynz][Y 1 ! [X 


Using node elimination method we retain only the buses 
where synchronous and induction machines are connected and 
eliminate the rest of the buses. For various cases, the compact 
models are developed in the following sub sections. 


[Zrr] = [Yrr] ! and 


4.1. SMES at wind park bus 


In D-Q component form, one can write Eq. (23) as: 


| Vino] = |GBipo| | Zip] + [GBipo] | Zipesm| (24) 
where: 


T 

[Ving] = [Vp] Vor Vp» Vaz Vps Vo] 
T 

[ogo] = [Ip Je: Ip2 Jaz Ip3 Io] 


[AposM] =[0 0 0 0 Ipsm Iosm]" 


Elements of GB,po are given in Appendix A. 

Eq. (24) can be combined with the stator equations of the syn- 
chronous and induction machines to obtain a compact model. For 
the sake of simplicity, all the induction machines can be repre- 
sented by a single equivalent induction machine run by a single 
turbine [9,10], on the assumption that all wind machines operate 
under same wind conditions. After transforming synchronous 
machine equations to D-Q axis, we obtain the compact model 
of Eq. (25). 


hpo = [Hipol ![Eipo — [G BipoMiposu] (25) 


where, the elements of [GB1po] and that of [H1po] are given in 
Appendix A and 


lno = Uni Tei Ip2 IQ2 Ips Ios]* 


T 
Emo = [Epi EQ E EQ Ep; EQ] 


liposM — [0 0 0 0 JIpsm Ios. 


4.2. SMES at load bus 


In D-Q (synchronously rotating reference frame) component 
form, for this case, one can write Eq. (23) as: 


[Vopo] = [GBzro| | pe] + [GBsro] | I2posw| (26) 
where 


T 
[Vno] = [Vp] Vor Vp? Vo» Vps Vol] 


T 
[bogo] = [Ip Jo: Ip2 Jaq2 lps Io] 


T 
[/oposM] = [/psM  JosM] 


Elements of GBzpg and [GB3pq] are given in Appendix A. 

Eq. (26) can be combined with the stator equations of the syn- 
chronous and induction machines as in previous case to obtain 
a compact model of Eq. (27). 


bpo = [Hopol [Exe — G BspolposM Q7) 
The elements of H2po are given in Appendix A. 
T 
[po] = [Ip Jo: !po lo? lps Io] 
T 
[E2pQ] = [Epi EQ Ep Et Ep; Eos] 


T 
[/2posM] = [/psM los] 


4.3. Wind park disconnection when SMES is at wind park 
bus 


In this case, we have, 


|Vspo] = |GBaro] [spe] |GBspo| | pes] Q8) 
where 


T 
[Vapo] = [Vp] Vai Vp» Vo2] 

T 
Lape] = [Ip Jo: Ip2 1o2] 


T 
[/3posM] = [/psM  JosM] 


Elements of GB4pQ and GBspo are given in Appendix A. 
The final compact model for this case is given by Eq. (29) 


pg = [Hno] [Espo — [G BapoMsposw] (29) 


The elements of [H3pqQ] and are given in Appendix A. 
4.4. Wind park disconnection when SMES is at load bus 


The model for this case is similar to the one obtained for 
Section 4.3. 


5. Simulation studies 


Based on the dynamic models of various components, pre- 
sented in the previous section, MATLAB programs were devel- 
oped and a range of simulations was performed. The behaviour 
of the system under various disturbances like wind park dis- 
connection, turbulent wind, sudden wind power variation and 
load changes is examined. Before running the dynamic programs 
for a particular disturbance, the load flow and initial condition 
program files are run to obtain steady state data prior to the 
disturbance (henceforth t< to). The power system considered 
comprises of two identical diesel units, each of 240 kW capacity. 
The wind park comprises of 10 identical wind turbines coupled 
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Fig. 4. (a) Frequency deviation (Hz) due to wind park disconnection, (b) voltage deviation at load bus (%) due to wind park disconnection, (c) active power supplied 


to SMES (W) and (d) reactive power supplied to SMES (VARS). 


to asynchronous generators, each of 12 kW. The load demand 
varies between 150 and 360 kW. The SMES unit considered 
requires a peak energy storage of just a few kJ and its converter 
rating is about 5% of the overall system capacity. 

The power quality results obtained with different configura- 
tions are compared for different cases. The SMES unit is first 
introduced at wind park bus and then shifted to load bus, to 
study the effect of the position of SMES unit, if any, on the 
power quality of the hybrid system. 


5.1. Casel: wind park disconnection 


The objective of this case study is to examine the behaviour 
of the power system with abrupt change in penetration level. 
Initially, the wind turbines are generating 100 kW and the power 
system load is 312 kW/200 KVAR, while the system frequency 
is 50 Hz. At t = to, the wind park along with associated capacitor 
banks is disconnected. According to Fig. 4(a and b), the power 
quality of the wind-diesel system is not acceptable [9]. 


For example the frequency deviation exceeds the acceptable 
limit of 1 Hz. The same figures reveal that with the introduction 
of SMES unit at wind park bus, the quality of the power supplied 
by the hybrid system is quite satisfactory. However, as revealed 
by the simulation results, power quality is not improved further 
as SMES unit is shifted from wind park bus to load bus. Fig. 4(c 
and d) shows active and reactive power supplied by SMES. The 
variations in these variables are due to continuous interaction of 
the SMES unit with the system during transients. 


5.2. Case 2: operation under turbulent wind 


The objective of this case study is to examine the behaviour 
of the hybrid power system when the wind turbines operate in 
turbulent wind. This mode of operation is very common for wind 
power based power systems. 

Initially, the wind turbines are supplying 50 kW of power and 
the load on the power system is taken as 260 kW/182 kVAR. For 
t» 19, the wind turbines operate in turbulent wind mode and 
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the power supplied by the wind turbines varies as shown in 
Fig. 5(a). 

Fig. 5(b and c) depicts the power quality results of the power 
system under consideration under turbulent wind operation. 
Although, for the wind-diesel system the maximum deviations 
in frequency and voltage are within limits, so is notthe case with 
rate of change of these variables [9]. This rate of change is an 
additional criterion for acceptable operation of power systems 
with wind energy conversion systems. For example, a high rate 
of change of voltages causes unacceptable flicker even if the 
maximum variation is in the neighbourhood of 2% [8,9]. Thus, 
the power quality ofthe wind-diesel system is not satisfactory for 
turbulent wind. The simulation results clearly establish the effect 
of introducing the SMES unit, on the power quality under most 
common mode of wind turbine operation. It can be inferred from 
the results that with the introduction of SMES unit, the maxi- 
mum deviations in system frequency and voltage are reduced to 
some extent. Although the power quality with SMES is satis- 
factory, the rate of change in load bus voltage is not reduced to 
a great extent because the fuzzy logic control proposed in this 
paper does not make use of rate of change of load bus voltage 
as one of the input variables. 


5.3. Case 3: operation under sudden wind power changes 


In this case study, the behaviour of the wind-diesel power 
system under step wind changes is studied. The wind power 
is considered to contain also a small stochastic component in 
the study. The stochastic component was generated by using 
MATLAB/SIMULINK toolbox and was then added to step wind 
changes thus resulting wind power profile of Fig. 6(a). 

Fort < to the wind turbines are supplying 50 kW of power and 
the load is 260 kW/182 KVAR. For t> fo the system operates in 
the wind power profile of Fig. 6(a). The simulation results for this 
mode of operation are shown in Fig. 6(b and c). These figures 
clearly show the positive impact of proposed control scheme 
on the power quality of the system, under sudden wind power 
changes. 


5.4. Case 4: load changes 


Change in load is a noteworthy source of perturbation in 
stand-alone power systems. This is because of the fact that a 
stand-alone power system has low inertia and is not able to con- 
tain the variations in system frequency and voltage unlike the 
utility grid. A change in static load is modeled as a change in 
network topology. 

Two different cases were studied as indicated in 
Tables 2 and 3. In one case, active power load changes are consid- 
ered while in second case the active as well as reactive power load 


Table 2 
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Fig. 5. (a) Wind power variation (kW), (b) frequency deviation (Hz) due to 
turbulent wind and (c) voltage deviation at load bus (%) due to turbulent wind. 
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Fig. 6. (a) Stochastic wind power variation (kW), (b) frequency deviation (Hz) 
due to stochastic wind disturbance and (c) voltage deviation at load bus (%) due 


to stochastic wind disturbance. 


Table 3 

Load profile for active/reactive power load disturbance 

Time (s) Active power (kW) Reactive power (kW) 
t<0 260 182 

0<t<10 360 200 


changes are considered. For t< to, the power supplied by the 
wind turbines is 50kW and the load on the system is 
260 kW/182 KVAR. In both the cases, wind power is considered 
constant. 

The simulation results pertaining to these two cases are 
depicted in Figs 7(a and b) and 8(a and b), respectively. It is 
observed from these results that the association of SMES unit 
with the wind-diesel system helps in maintaining the power qual- 
ity in this case also. 

In none of the case studies, the SMES coil current exceeds 
500A. Thus an SMES unit with a peak-stored energy of just 


Frequency Deviation (Hz) 
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a 
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Fig. 7. (a) Frequency deviation (Hz) due to active power load disturbance and 
(b) voltage deviation at load bus (%) due to active power load disturbance. 
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Fig. 8. (a) Frequency deviation (Hz) due to active and reactive power load dis- 
turbance and (b) voltage deviation (96) due to active and reactive power load 
disturbance. 


a few kJ is sufficient to deal with all these disturbances. The 
examination of active and reactive power curves of the SMES 
unit reveal that the SMES converter rating required is just 596 
of the overall system capacity. 


6. Conclusions 


An effective control scheme, based on SMES unit, for power 
quality improvement of a typical multimachine wind-diesel 


power system is proposed. The effect of insertion of SMES 
unit at two alternative places is explored. Complete models 
for various configurations and case studies are developed. The 
SMES unit is operated in four quadrants using a fuzzy logic 
control strategy. The proposed scheme is investigated for its 
impact on the system performance for various disturbances like 
wind park disconnection, wind gusts, abrupt wind changes and 
load perturbations. The simulation results show that the pro- 
posed scheme helps to maintain the power quality of the sys- 
tem, within the limits, for various perturbations. The simulation 
results also reveal that the quality of power is not improved 
much as the SMES unit is shifted from wind-park bus to the load 
bus. 


Appendix A 
A.1. System data 
Base power 400 kVA for the whole system. 
A.2. Synchronous generators 
Nominal power (kW) = 240, nominal voltage = 380 V. 
Nominal apparent power (kVA) = 240, rs = 0.067, Xa = 260%, 
X42156, X4-2296, X4-1596, P de =22.5%, T4,-1.182s, 
T =0.015s, Too =1.18s. 
A.3. Induction generators 
Nominal power (kW) = 12, nominal voltage = 380 V, nominal 


speed = 1545 rpm, rs = 0.022 p.u., xs 20.181 p.u, xm = 1.605 p.u., 
r; = 0.022 p.u., x, 2 0.181 p.u., HA =0.259 s. 


A.4. SMES 
Lsm = 0.264 H, 19, 2350 A, V9, 2 42.5|V|, Ke 0.01 
Elements of GBjpo are given by 
[GBipo] 
Gli —Bly Gli —Bli Gly —Bli3 
Bly Gly Bl2 Gly Bly Gli 
Glo, —Bla Glo, —Blo Glo —B1x% 
= Bla Gly Blo Gla Bliz G123 
Gl3; —Bl3; Gl3. —Bl32 Gla3 —B133 
Bl3, G13 B132 Gl32 Bl3}  Gla3 
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HipQ(4. 3) = X(n cos? à» + Xt sin? ô2 + B122 
Hipo(4, 4) = rj + (Xp — X) sin ó? cos 62 + G125 


HipQG, 5) = r3 + G133 
Hipo(, 6) = —(X3 + B133) 


Hipo(6, 5) = X} + Blas 


Elements of [Hipo] 
Hipo(l, 1) = rs + (Xi — XG sin 61 cosó; 4 Gly  Hipo(4, 1) = Bla 
Hipo(1,2) = — (X7, sin? 8, + X% cos?ài-- Bly) — Aipg(4,2) = Gla 
Hipo(1, 3) 2 Gli? 
Hipo(1, 4) = —Blio 
Hipo(1, 5) = Glis Hipo(4, 5) = B123 
Hipo(1, 6) = — B113 Hipo(4, 6) = G123 
HipoQ. 1) = X7, cos? 3, + X"; sin? 3j + B111 Hipo(5, 1) = G31 
Aipg(2, 2) = rs + OG a X'i) sin 61 cosó; 4 Gly,  HipQ(5, 2) = — B131 
Aipg(2, 3) = Bli? AipQ(5, 3) = G132 
Hippol, 4) = Gli? HipQG, 4) = — B132 
HipqQ. 5) = B113 
Aipg(2, 6) = Gl13 
Hipo, 1) = Glz; Hipo(6, 1) = B131 
Hipo, 2) = — Bla) Hipo(6, 2) = G13, 
Ai pg(3, 3) = rj + (Xo — X) sinó? cosd2 + Glo? HijpQ(6, 3) = Bla; 
Hypa. 4) = — (X7, sin? 35 + Xj, cos? 62+ Bln) ^ HipQ(6.4) = G32 
HipoG, 5) = Gl73 
HipoQQ, 6) = — B123 


The elements of GBzpg and GBspo are given by 


G211 
B211 
G221 
B221 
G231 
B231 


[GBopo] = 


Hipo(6, 6) = rs3 + G133 


-B211 G22 —B2ı2 G2: — B213 G214 B2i4 G2əą_  B2344  G234 B234 i 
G2j; B22 G22 B23  G215 [GBano] = —B24 G2\4 -—B24 G24 —B234 G234 
— B21 G2» —B2n G23 -—B2»5 

G21 B2» | G222 B23 G23 

— B23, G232 —B232 G233 —B233 

G21 B2» G2» B2: G23 Elements of [Hpo] 


Hopo(1, 1) = rs1 + (XG, — Xt) sin ôi cosó1 + G2j1 
H»po(1,2) = — (OX sin? 81 + X4] cos? 81 + B211) 
Hopo(1, 3) = G212 

Hopo(1, 4) = — B212 

Hopo(1, 5) = G213 

Hopo(1, 6) = — B213 

H»poQ. 1) = Xs cos? ô; + XG, sin? 61 + B211 
Aapg(2, 2) = ra + OG — X) sin; cosó; + G21 
HopQQ. 3) = B212 

HopoQ., 4) = G212 

HopQQ2. 5) = B213 

HopoQ., 6) = G213 

HopoG, 1) = G221 

Apg(3, 2) = — B22; 

Hop9 (3, 3) = rs2 + (XG — Xp) sin 52 cos ô2 + G22 
Hagl, 4) = — Qt sin? 8; + X/ cos? 85 + B222) 
Hpg, 5) = G223 

H5poQG, 6) = — B223 


Hopo(4, 1) = B221 

Hopo(4, 2) = G22) 

Hopo(4, 3) = X? cos? 6 + X sin? 59 + B225 
Hopo(4, 4) = rs2 + (Xp — X'p) sin ô2 cos 62 + G222 
Hopo(4. 5) = B223 

Hopo(4, 6) = G223 

HopoG, 1) = G231 

Hopo(G, 2) = — B231 

FApQ(5, 3) = G232 

ApQ(5, 4) = — B232 

Hopo(, 5) = rs3 + G233 

Hopo(5, 6) = —(X5 + B233) 

Hopo(6, 1) = B231 

Hpo(6, 2) = G23; 

Anna (6, 3) = B232 

Hopo(6, 4) = G232 

Hopo(6, 5) = X} + B233 

Hopo(6, 6) = rs3 + G233 
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The elements of [GB4po] and [GBspo] are given by 


G31; —B31; G32 —B312 

B311 G31; B32  G3p 
[GB4po] = 

G32 —B321 G32 —B322 

B32, G32 B3» G3» 

G313 B33 G33 B3» 
[GBspo] = 

— B313 G33 —B3»; G35 


Elements of [H3po] are given by 
Hspo(l, 1) = rsi + (XG, — Xt) sin ôi cos 6, + G311 
Hspo(1, 2) = — (X7; sin? 81 + X/, cos281 + B311) 
Hapo(1, 3) = G312 
Hapo(1, 4) = — B312 
HspoQ. 1) = X, cos? 8) + X sin? 5) + B311 
A3pQ(2, 2) = rsı + (Xo — 
FA3pQ(2, 3) = B312 
FA3pQ(2, 4) = G312 


4p sin; cos 6; + G311 
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